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-L: Introduction. Characteristics of suprathermal particles accelerated 

by quasi-parallel interplanetary traveling shocks have been generally 
explained in terms of a first order Fermi mechanism (e.g. see review in 
ref. 7). Such models require diffusive scattering of particles upstream 
of the shock. This scattering is characterized by a local diffusion 
coefficient, k, which is determined by the local power density of waves 
in the upstream region as described by Lee (7). A number of studies 
have investigated the behavior of k of £1 MeV upstream ions close to 
interplanetary shocks. Scholer et al. (8) have used results of first 
order Fermi shock acceleration theory (e.g. ref. 1) to derive the 
diffusion coefficient and its energy dependence from the measured 
gradient of the upstream particle intensity. Van Nes et al. (9), using 
the same approach, obtained the spatial and energy dependence of k for 3 
interplanetary shocks over a more extended energy and distance range. 
Klecker et al. (6) deduced the spatial dependence of k in the upstream 
region for 30 and 130 keV protons from the measured first-order 
anisotropies and intensity gradients. 

In this paper we examine the dependence of the diffusion coef- 
ficient of suprathermal upstream protons on distance from the November 
12, 1978 interplanetary traveling shock using a different approach. 

Unlike previous studies our method, which is based on measurements of 
particle streaming and intensity gradients, does not rely on predictions 
of shock acceleration theories or require first-order expansions. We 
have chosen to examine the local spatial variations of k upstream of the 
November 12, 1978 shock because the characteristics of this quasi-paral- 
lel shock have been extensively studied (e.g. ref. 5), and also because 
of its favorable geometry (i.e. B field nearly radial). The initial 
results of this study have been reported by Gloeckler et al. (3). 

2^1 Instrumentation and Method of Analysis. For this study we use the 

counting rate data from the University of Maryland/Max-Planck-Institut- 
Garching ULECA sensor (see ref. 4 for details) on ISEE-3 which was 
placed in a halo orbit at a radial distance of ~230 Rg upstream of the 
earth. Of relevance to the present discussion are the capabilities of 
the electrostatic deflection vs^ energy ULECA sensor to reliably separate 
protons from alpha particles and to determine their differential inten- 
sities in three energy bands centered on energies E of 33, 66 and 132 
keV/e, and respective band widths AE/I of 20, 27 and 37%. The sensor 
has a 60 fan-like acceptance aperature with its wide angle in a plane 
perpendicular to the ecliptic. The three proton counting rates are 
available every 128 sec in each of eight 45^ sectors of the ecliptic 
plane. 


For the present analysis we use the 33 and 66 keV sectored proton 
counting rate data to derive the proton distribution function, f, in the 
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the two segments of f(v) join 
smooth, as shown in Figure 1. 
requires no prior knowledge o: 
spectrum nor need V F /v be small. 


Figure 1 . Rest frame proton phase 
space density, f(v). For V F = 999 
km/s and ^ = 177° segments of f(v) 
from the 33 keV proton rate channel 
(open squares) and the 66 keV 
channel (solid circles) join 
smoothly to form a common spectrum. 


frame of reference (rest frame) in 
which f is assumed to be isotropic, 
and find the velocity, \T F » of that 
rest frame relative to the space- 
craft. We follow the technique of 
Gloeckler et al. (2) for mapping 
the sectored counting rate data 
into segments of the rest frame 
distribution function f(v), where v 
is the particle speed in the rest 
frame, making use of the invariance 
properties of f • When the proper 
rest frame velocity, ^ F , is chosen, 
up and the resulting spectrum becomes 
We note in particular that our method 
the shape of the rest frame energy 


3. Observations. Using a least-squares technique we derive for^each 
available 128 sec interval the rest frame velocity and the rest frame 
distribution function (over a range of particle speeds from -1000 to 
5000 km/s). Figure 1 shows an example of a 128 sec averaged distri- 
bution function f(v) obtained ~ 4-10 10 cm upstream of the shock. We 
note that the shape of f(v) tends to be between that of a pure power law 
and an exponential, and that the spectrum becomes harder with distance 
from the shock (3). In the present analysis an exponential form is 
assumed and the e-folding speed, v Q , is computed from a linear fit to 
the data as shown in the figure. 


In Figure 2 we plot the derived speed, V F , (solid circles) and 
direction of motion, (upper panel) of the suprathermal proton rest 

frame for each available 128 sec interval as a function of time as the 
shock is approached. Notice that the direction o F & 

within ~10° antisunward. Also shown is the sum of the solar wind, Vg„, 
and Alfvdn, V., speeds taken (or derived) from data published in 5). 
While far ahetd of the shock V p exceeds (V gw + V A ) by as much as 1000 
km/s, close to shock this difference is about a few hundred km/s, and 
behind the shock the two velocities are the same within expenmenta 
errors, implying pure convention at a speed equal to the sum of solar 
wind and Alfvfen speeds. If we interpret the difference between the 
speeds V v and (V sw + V A ) upstream of the shock to be the diffusive 
streaming of suprathermal ions along the magnetic field B (near y ra la 
for this shock) which results from the observed upstream particle 
density gradient, the local diffusion coefficient along B, k„ (^.v) at 
a distance r. = (t„ u -t< )V aU from the shock (V gh and t gh are the shock 
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speed and arrival times respectively) may be related to quantities 
measured at r^ as follows: 

(1) £ = 0 = fy v (-||) ]*e B • [^ f (v)-(^ sw +\T a ) ] + K || (r i ,v)*e B .^f(v). 

$ is the differential streaming (which in zero is the particle rest 
frame), and is the unit vector along Since for the November 12, 
1978 shock $ and were nearly radial, and assuming that was 

directed along $ in the upstream region, and that f(v) has an exponen- 
tial dependence on v with e-folding speed v Q , eq. (1) may be simplified 
to 



To obtain < at a distance r^ upstream of the shock values for v are 
obtained from the slope of the ith distribution function, and° the 
particle density gradient is determined from the neighboring values of f 
at v - 3100 km/s (50 keV). V p and (V A + V gw ) are those plotted in 
Figure 1 at times t ± related to r. by r ± «( t s h“ t i) V sh 5 where we used 
values for V gh of 612 km/s and t gh of 0:28:18 as given in (5). 

Figure 2 , Speed and 

direction of motion of 
the rest frame (solid 
symbols) and the solar 
wind plus Alfvfen speed 
(V gw + V A ) vs^ time be- 
fore and immediately 
after shock passage. The 
difference between Vp and 
(V gw + V A ) is interpreted 
in terms of diffusive 
streaming. Behind the 

shock V F ~(V SW + V A ) 
implying pure convection 
of 50 keV protons at that 
speed • 

4. Spatial Dependence of 

Kjj » In Figure 3 we show 

values of K (j for ~ 50 keV 

protons as a function of 

distance from the shock, 
where k [( is computed 

using eq. (2). The un- 

certainties shown come 
primarily from inaccuracies in determining the local density gradients. 
We note that there may be systematic errors (£30%) in the values of k 
resulting from using the simpler eq. (2) rather than eq. (1) which in- 
corporates vector quantities. Consistent with expectations (7) we find 
that k increases exponentially from a value of -7 *10 17 cm 2 /s near the 
shock to -2.5 *10 18 cm 2 /s at - 4.5 *10 4 * * * * * 10 cm from the shock with an e- 
folding distance of 3.4«10 18 cm. What is surprising is the subsequent 
decrease (again exponential, with e-folding distance of -1.1 *10 10 cm) in 
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the value of k with distance from the shock. 

5. Di scussion and Conclusions . The factor of ~ 4 increase in k over a 
distance of ~ 4.5»10 iU cm is consistent with a comparable decrease in 

the^ower spectral density of -ayes at I-IO'* H. («-««.» of ref. . 

The value of 7*10 i0 cm^/s of 
keV) near the shock is in excellent 
agreement with prediction (for this 
shock) based on a self-consistent 
theory for wave excitation and particle 
acceleration upstream of interplanetary 
shocks (7). Our values of k and its 
spatial dependence are also in quali- 
tative agreement with those reported in 
(6,9). What is puzzling is the de- 
crease of k (50 keV) at larger 
distances which is not predicted by 
current theory. To pursue this point 
further we computed < ^ (100 keV) fol- 

lowing the procedure outlined above but 
now using the 66 and 132 keV proton 
counting rate data. We found that the 
rest frame speeds V-p characteristic of 
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Figure 3. Spatial dependence 
of k on upstream distance 
from Ihe shock. 
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~100 keV protons are smaller than those 
shown in Figure 1 for 50 keV protons. 
Ratios of K100/K50 were determined to be ~3.0, 1.2 and 1.8 near the 
shock at r = 4 *10 10 cm and 8*10 10 cm respectively. Assuming a simple 
power * law dependence of k on particle energy these values correspond to 
*1 a e +1»6 e^ 0 * 26 , E +0 * 85 respectively. We therefore conclude in agree- 

ment with’ (9) that one cannot characterize the energy dependence (or the 
spatial dependence) of K in a simple way. Oar results indicate an in- 
crease in the wave power density beyond r -5-10 cm and in general a 
more effective scattering of 50 keV protons compared to 100 keV 
protons. This hypothesis can be checked using local wave power spectrum 
measurements as a function of distance from the November 12, 1978 
interplanetary traveling shock. 
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